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We present density-matrix renormalization group results 
for the ground state properties of two-leg Hubbard ladders. 
The half-filled Ifubbard ladder is an insulating spin-gapped 
system, exhibiting a crossover from a spin-liquid to a band- 
insulator as a function of the interchain hopping matrix ele- 
ment. When the system is doped, there is a parameter range 
in which the spin gap remains. In this phase, the doped holes 
form singlet pairs and the pair-field and the "4^^" density 
correlations associated with pair density fluctuations decay as 
power laws, while the "2^^" charge density wave correlations 
decay exponentially. We discuss the behavior of the expo- 
nents of the pairing and density correlations within this spin 
gapped phase. Additional one-band Luttinger liquid phases 
which occur in the large interband hopping regime are also 
discussed. 



I. INTRODUCTION 

There are a number of new materials which contain 
weakly coupled arrays of metal-pOxide ladders. For ex- 
ample, SrCu203El and La2Cu2O5&0 contain two-leg lad- 
ders with Cu-O-Cu rungs, and (V0)2P207O contains iso- 
lated two-leg V-0 ladders. Alternatively, Sr2Cu305 con- 
tains three-leg ladders and Sr„Cu„+i02n+i has n-|-l-leg 
ladderaJl. Magnetic susceptibility ancj-|nuclear relaxation 
time measurements provide evidenceO that in the insu- 
lating state, the even-leg ladders have a spin gap while 
the odd-leg ladders have gapless spin excitations. It is 
of great interest to understand what happens when holes 
are doped into these systems, and, in particular, what 
happens when holes are doped into the spin-gapped two- 
leg ladders. Recently, Hiroi and TakanoH reported that 
the two-leg ladder compound La2Cu205 could be hole- 
doped by substituting Sr for some of the La. They found 
that the conductivity of the doped ladders has a metallic 
behavior and that magnetic susceptibility data suggests 
that a spin gap remains in the lightly doped system. 

Here we examine the properties of a Hubbard model 
on a two-chain ladder in order to understand the ground 
state-behavior of the undoped and doped two-leg ladders 
in terms of an itinerant electron model. We will show 
primarily numerical results for the energies and equal- 
time correlations of the ground and low-lying states ob- 



tained usiHg the Density Matrix Renormalization Group 
techniqueLJ (DMRG). -While parts of this work have been 
published previouslyQ, here we will discuss the results 
in more detail, make comparisons to weak- and strong- 
coupling analytic pictures, and also present some new 
results on the local structure of the pairing, the ^%kp^^ 
charge-density-wave (CDW), and on the behavior of the 
spin and charge gaps in the quarter-filled system. 

The Hubbard and related models on a two-leg lad- 
der have received much theoretical attention recently, via 
various analytic approximations, as well as a variety of 
numerkial techniques. We direct the reader to a recent 
reviewcl and the references contained therein. Previous 
numerical work includes-ianczos calculations for a t-J 
model by Dagotto et al.ll2l who suggested that antiferro- 
magnetic 5 = 1/2 coupled chains should have a spin gap 
and that if the two-chain system could be doped, it would 
have enhanced superconducting or charge-density-wave 
correlations in the ground state. In Ref. 11, the Heisen- 
berg model was treated using exact diagonalization and 
analytic techniques in order to determine the ground 
state properties, the temperature dependence of the spin 
susceptibility, . a i adi the spin excitation spectrum. Rice 
and coworkeral^llj discussed a t-J model and suggested 
that if this system were lightly doped, the spin gap phase 
would persist and a ground state with dominant super- 
conducting correlations could be realized. The one- and 
two-particle dynamical correlation functions were calcu- 
lated for the half-filled Hubbard ladder using Quantum 
Monte Carlo and a Maximum Entropy analytic continua- 
tion in Ref. 14. There have also been a number pf weak- 
coupling renormalization group calculationaijE3 which 
find evidence for a variety of phases. We will discuss the 
results of one of these calculationscJ in more detail below 
and compare them to our numerical calculations. 

In the following, we present numerical evidence that 
indicates that the half-filled two-leg Hubbard ladder is 
a spin-gapped insulator for all non-zero values of the 
Coulomb interaction U, and the interchain coupling t^. 
We show that for weak [/, there is a well-defined spin- 
liquid insulator to band insulator transition that corre- 
sponds to the [/ = one-band to two-band transition. 
For stronger [/, this transition evolves into a smooth 
crossover. When the system is doped, the charge gap dis- 
appears, but the spin gap remains for a range of t±^ asso- 
ciated with overlapping bonding and antibonding bands. 
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For t± larger than this range, the numerical results are 
consistent with Luttinger liquid behavior. Within the 
spin gap phase, we show that there is an attractive pair- 
ing interaction, a d-wave-like structure of the pair wave 
function, and algebraically decaying pairing and "4fci?" 
CDW correlations. However, the CDW correlations do 
not seem to decay in the manner predicted by., weak- 
coupling bosonization treatments of the systemEilo. 

The organization of the paper is as follows: in sec- 
tion H, we review the basic properties of the model and 
discuss the predictions of weak-coupling renormaliz ation 
group and bosonization calculations. In section [I A, 



we exhibit numerical results for the half-filled system, 
discuss the properties of the spin liquid state and the 
crossover to a band insulator as a function of the perpen- 
dicular hopping. In section II B, we explore the effects 



of doping holes into the system, and discuss the phases 
present as a function of the perpendicular hopping at 
various fillings. We discuss in detail the behavior of the 
holes within the gapped spin liquid state that persists as 
the isotropic system is doped. 



II. THE TWO-CHAIN HUBBARD MODEL 

The single band Hubbard model on two coupled chains 
of length L has the Hamiltonian 



2, Act 



h.c.) - ti_ ^(4i^c. 2,, + h.c.) 



(1) 



Here c- and create and destroy, respectively, an 
electron on rung i and chain A with spin a, and xa — 



c\ )^ Ci xcr- We will also use the notation rii 



The hopping integral parallel to the chains is t, the hop- 
ping between the chains t± , and U is the on-site Coulomb 
interaction. For the remainder of this work, we set t = 1 
and measure all energies in units of t. Within this pa- 
per, we will discuss a system with open boundary condi- 
tions both between the two chains and at the ends of the 
chains. In other words, site i, 1 is connected through only 
one hopping term to site z, 2, and site 1, A is not directly 
connected to site L, A. 

The non-interacting, [7 = 0, Hamiltonian can be 
diagonalized by writing the hopping term in terms of 
bonding and antibonding states on a rung and fourier- 
transforming parallel to the chains. The energy for the 
infinite two-chain system is then given by 



— — (2t cos k + t± cos fc ) 



(2) 



with k = (k,k±), where fcj_ = and fcj_ = tt corresponds 
to the energy of the bonding and antibonding band, re- 
spectively, and k is the momentum along the chains. For 
t± = 0, the two bands will be degenerate, and for nonzero 
t±, they will be separated by an energy 2t±. 



Due to the band structure, the = phase diagram 
in the t± - (n) plane, shown in Fig. 0, exhibits significant 
structure. Both bands will be occupied when tj_ < tj_c 
(shaded region), whereas when t± > t±c, only the bond- 
ing band will be occupied. Here tj_c = 1 — cos7r(n) and 

is shown by a solid line and (n) = (X^o- '^1 x a'^i,^-'^) 
band filling. At half-filling t±c — 2, and the system is a 
two-band metal for t± < 2 and a band insulator with a 
completely occupied bonding band for t± > 2. 




<n> 

FIG. 1. The (n) - t± phase diagram for the L'^ = system. 
In the shaded region, both the bonding (lower) and antibond- 
ing (upper) bands are occupied. At the solid line, the lowest 
part of the antibonding band just touches the Fermi surface, 
and on the dashed line, the bonding band is half-filled. 

Along the dashed line in Fig. |l|, the bonding band is 
half-filled. At half-filling {(n) = 1), this occurs at t± = 
0, where both bands are half-filled. As (n) is reduced, 
the t± at which the bonding band is half-filled becomes 
larger and the occupation of the antibonding band at this 
point becomes smaller. When (n) — 0.5, the bonding 
band is half-filled only when the antibonding band is 
completely unoccupied, i.e. for t± > t±c- 

Starting from the band-structure of the noninteract- 
ing system, one can treat the system within a weak- 
coupling picture by linearizing the band structure around 
the Fermi points and taking the continuum limit. The 
resulting model can be treated within a renormalization 
group j-oi- within a bosonization picture. Balents and 
Fisheir3 treated the small U, but arbitrary t±_ limit in a 
systematic, controlled procedure. In this approach, weak 
coupling RG equations were numerically integrated for 
each value of (n) and t± , evolving from infinitesimal U to 
small, but finite U. At this point bosonization methods 
were used to analyze the resulting Hamiltonians. Their 
weak-coupling results for the phase diagram agree sur- 
prisingly well with our results at rather strong coupling. 
For this reason we will adopt their notation and compare 



2 



our results with the results of their approach, which we 
will denote BFRG. An alternative strong-coupling ap- 
proach, which is less predictive but perhaps more intu- 
itive, is based on a s h ort-range resonating-valence-bond 
variational ansatzC3E3. This approach works particularly 
well in explaining the structure of the pair wavefunction 
in the doped, spin-gapped phase. We will discuss this in 
more detail in Section II B. 



Since, in a weak-coupling picture such as that used by 
Balents and Fisher, there are in the general case four 
Fermi points (at ±kp and ±fc^), a bosonization treat- 
ment yields four possible modes: two spin modes, one 
symmetric and one antisymmetric with respect to inter- 
change of the chains, and two charge modes, one symmet- 
ric and one antisymmetric. This is a generalization of the 
bosonization picture for the one-dimensional Luttinger 
model, which results in one spin and one charge mode. 
In general, each of the four modes can either be massive 
(i.e. gapped), or massless (with gapless excitations). Ba- 
lents and Fisher classify the possible phases according to 
the number of massless charge modes ric and spin modes 
Us, using a notation CricSns- In the two-chain system, 
there are nine possible phases, ranging from a COSO phase 
with excitation gaps in all four modes to a C2S2 phase, 
with all four modes gapless. For example, the two-chain 
analog of a Fermi liquid is a C2S2 phase, which occurs at 
[7 = when both bands are partially filled. The BFRG 
calculation, which is valid for weak U but arbitrary tj_, 
finds seven different ground-state phases in the {n)-t± 
phase diagram. 

Since the origin of the BFRG phases can be under- 
stood qualitatively in the context of the U = phase 
diagram, we will briefly describe them here. At half- 
filling, the system is in a COSO (i.e. spin and charge 
gapped) phase for all [/ > and t± > 0. For t± < 2.0, 
umklapp processes involving two particles scattering be- 
tween the Fermi points lead to gaps opening in all four 
modes, producing a spin liquid insulator. For t± > 2.0, 
the system is a band insulator, as in the U — case. 
Upon hole doping for t± > t±c, only the bonding band 
comes into play and the behavior is that of a one-band 
Luttinger liquid. Within this phase, there are gapless 
spin and charge modes, i.e. a ClSl phase, except at 
quarter-filling, where the relevant bonding band is half- 
filled and umklapp processes within the band lead to a 
charge-gapped COSl phase. For t± < t±c and (n) < 1, 
the BFRG yields a CI SO phase for most of the phase 
diagram. This phase, a spin liquid with one gapless sym- 
metric charge mode, was also predicted by other, earlier 
weak-coupling renormalization gr oup i CaJ culationsE^IE^, 
by various strong-coupling pictureal3'E3E3, and has been 
found withip .thfl-feibbard and t-J models by numerical 



calculationaiilliij'EJEj. Other phases are found in certain 
regions where special processes become relevant. One 
such region is in the vicinity of t± = t±c^ where the bot- 
tom of the antibonding band just touches the Fermi sur- 
face so that its dispersion must be treated quadratically, 
rather than linearized. Here the BFRG calculation finds 



a CISO phase on the band-transition line. For t± slightly 
below t±c, the Fermi velocities in the bands are quite dif- 
ferent and the BFRG finds narrow regions of first C2S2, 
then C2S1 phases as t± is reduced. Umklapp processes 
become relevant where the bonding band is half-filled, 
along the dashed line in Fig. |^. Here the BFRG calcu- 
lation finds a C1S2 phase along this line near half-filling 
and near quarter-filling, with an intermediate region in 
which the CISO phase remains present. 

It is interesting at this point to compare the behavior 
of the two-chain system to that of the one-dimensional 
interacting electron gas (the ID Hubbard model, for ex- 
ample). For repulsive short-range-jnteractions, the ID 
electron gas is a Luttinger liquidcJ. A Luttinger liq- 
uid has gapless spin and charge excitations (a ClSl 
phase), except when umklapp processes become relevant. 
The dominant long-distance behavior of the spin-spin, 
density-density, and on-site s-wave correlation functions 
is power-law decay, and to leading order isllj 



(S(r) ■ S(0))tl = r-(i+^^) cos(2/ci.r) 
(n(r)n(0))TL = r"(i+^''' cos(2fcFr) 



(A,(r)At(0))TL = 



(3) 



where S(r) is the total spin at site r, n{r) is the electron 
density, and Aj(r) creates an on-site s-wave pair. The 
exponent Kp is a non-universal parameter dependent on 
the system parameters. There are logarithmic correc- 
tions multiplying these expressions which have been left 
out for simplicity. Umklapp processes become relevant 
at half-filling for the Hubbard model, at which point 
a charge gap develops, leading to a COSl mode. The 
strong-coupling limit of the half-filled Hubbard model 
is the Heisenberg model which has no charge degrees of 
freedom, and in ID is known to have gapless spin excita- 
tion^. The spin-spin correlation decays to leading order 



(S(r) • S(0))Hcis = cos(^r). 



(4) 



In the half-filled Hubbard model, the spin-spin correla- 
tion has this form and other correlation functions decay 
exponentially. According to the BFRG calculation, the 
doped two-chain system should behave as a Luttinger liq- 
uid when t± > t±c, and the filling of the Luttinger liquid 
corresponds to the filling of the bonding band. Therefore, 
the correlation functions should decay as power laws gov- 
erned by one exponent, as in Eq. (|^) in this regime, and 
by Eq. (^ at quarter-filling. 

For attractive short-range interactions, the ID elec- 
tron gas falls into the Luther-Emery universality classES. 
Within the Luther-Emery model, there is a spin gap and 
one low-lying charge mode (CISO). The s-wave pairing 
correlation function and the CDW correlation function 
both decay algebraically, with a power parameterized by 
the CDW exponent Kp. The leading behavior of these 
two correlation function is given by 
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(n(r)n(0))LE = r'^" cos(2fcFr) 
(A,(r)At(0))LE-r-i/^^. (5) 

Here Kp is also dependent on the model parameters, such 
as the band-filling, or the strength of the Coulomb inter- 
action, \U\, in the attractive Hubbard model. Within 
the BFRG calculation, the doped two-leg ladder falls 
into a Luther-Emery-like CI SO phase for much of the 
region for t±_ < t±c- It is therefore a relevant question 
whether there are correlations analogous to those in Eq. 

with a reciprocal relation between the exponents in 
this phase. Roapaization treatments of the two-chain 
system predictE2l'E2l that the relevant CDW operator will 
occur at the "4/ci?" wave vector q* — {2[kp + kp], 0) with 
an effective density operator = ni^ini^2, and that the 
pair wave function will no longer be s-wave, but c?-wave- 
like with the largest amplitude across a rung. Since the 
behavior of the correlation functions is dependent on the 
parameter Kp, which is a non-universal quantity which 
can depend on the parameters of the model, it is im- 
portant to determine the exponents within a particular 
phase. The least restrictive criterium for the supercon- 
ducting pairing being dominant over the CDW is that 
the correlation function decays more slowly at long dis- 
tances, which occurs when Kp > 1. However, the consid- 
eration of the effects-of impurity scattering and crossover 
to three dimensionscJ can lead one to the more stringent 
requirement that Kp > 3, in order to obtain a stable 
superconducting state. 

The numerical results shown in this work are all cal- 
culated using the DMRG techniqueQ on finite lattices of 
2 X 8 to 2 X 64 sites. We obtain the energies and equal- 
time correlation functions of the ground state and the 
low-lying excited states of the finite cluster. While the 
DMRG technique gives energies that are, in principle, 
variational, it has proven to give quite accurate results 
for ID quantum lattice systems. The method provides a 
controlled way of numerically diagonalizing a finite sys- 
tem within a truncated Hilbert space, by successively 
building up part of the system using a real-space block- 
ing transformation, and then using the reduced density 
matrix to truncate the Hilbert space of that part of the 
system in a controlled way. One can increase the ac- 
curacy by increasing the number of states kept, and can 
examine the convergence with the number of states. Here 
we typically keep 400 states per block, although in the 
numerically more difficult cases, such as the calculation 
of the spin gap on doped 2x32 lattices with small t± , we 
keep up to 550 states. Truncation errors, given by the 
sum of the density matrix eigenvalues of the discarded 
states, vary from 3.7 x 10~^ in the worse case to O(10~®) 
in the best cases. This discarded density matrix weight is 
directly correlated with the absolute error in the energyQ. 
It was only possible to obtain accurate calculations on the 
larger lattice sizes, 2 x 40 and 2 x 64, in certain parame- 
ter regimes, in which the convergence with the number of 
states was relatively rapid. In other regimes, the largest 
lattice size was 2 x 32. We estimate that the maximum 



errors on the quantities shown in this paper are at most a 
few percent and typically are of the order of the plotting 
symbol size or less. We apply open boundary conditions 
to the lattice because the DMRG method is most ac- 
curate for a given amount of computational effort with 
these boundary conditions. 

A. Half-filling 

As discussed above, the undoped two-leg ladder ma- 
terials are gapped spin-liquid insulators. Analytic and 
numerical studies of the Heisenberg model, which is the 
strong-coupling limit of the Hubbard model at half- 
filling, determine that th e sys item is a spin-liquid insula- 
tor for isotropic couplingalJO'EHl. There is also numerical 
evidence that this spin-liquid state is pr|eseii|t-£pr all ratios 
of the interchain to intrachain couplingi2l'E3E3. However, 
it is not clear how the crossover from the weak-coupling 
picture dominated by the band structure effects to the 
strong-coupling spin-liquid state takes place within the 
Hubbard model. We will explore this question in this 
section. 

In order to see the structure of the correlations at half- 
filling, we first examine the spin-spin correlation function 
defined as 

S{z,j,X)^{Ml^Ml,) (6) 

where M^^ = ?^i,At ~ '^i.AX is the z-component of the on- 
site magnetization. This correlation function, multiplied 
by a factor (—1)' which removes the antiferromagnetic 
g = TT structure, is shown in Fig. ^ on a semilog scale. 
Since the correlation functions are straight lines on the 
semilog plot, the decay is exponential, as we would ex- 
pect in a gapped spin liquid. The spin-spin correlation 
function decays more rapidly as U is increased, implying 
the correlation length decreases with U. The behavior of 
the correlation length, obtained by fitting a straight line 
to the curves between 6 and 19 lattice spacings, is shown 
in the inset. The dashed line in the inset is the spin 
correlation length for the two-chain isotropic Heisenberg 
system, ^ = 3.19, taken from Ref. 23. As can be seen, 
the correlation length for the Hubbard model converges 
to the Heisenberg value for large U, as one would expect. 

We have examined the spin-spin correlation function 
for a variety of U and t± at half-filling, down to t± = 0.2 
at [/ = 1, and f/ = 0.5 at = 1.0. The decay of the 
correlation functions is consistent with exponential decay 
and a finite correlation length at half-filling for all t± 
and U we have examined. Therefore, we believe the spin 
liquid phase is present at all t± and U at half-filling. This 
is consistent with weak-coupling renormalization group 
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treatmentslliH. 




FIG. 2. The spin-spin correlation func- 

tion {-lY (Ml^Ml-^) for various U at half-filling ((n) = 1.0) 
plotted on a semilog scale. The inset shows the correlation 
length extracted from the slope of the lines as a function of 
U, with the dashed line the Heisenberg value from Ref. 23. 

We have also examined other correlation functions, 
such as the density-density and pair correlation func- 
tions, for the half-filled system. In all cases they also de- 
cay exponentially, with a correlation length shorter than 
that of the spin-spin correlation function. 

We next consider the charge and spin gaps defined by 

Ac = [Eo{N - 1, iV - 1) + Eo{N + 1,N + 1) 

~2Eo{N,N)]/2 (7) 

and 

As = EoiN + l,N-l)-EoiN,N), (8) 

respectively. Here Eq{N^ , N^) is the ground state en- 
ergy for iVf spin-up electrons and Ni spin-down elec- 
trons. Since our current DMRG program cannot target 
states of a particular transverse symmetry, we can only 
calculate the gaps to the lowest-lying spin and charge 
excited states, irrespective of symmetry. We show spin 
and charge gaps for a range of U and t± in Fig. ^. For 
U ~ 0, the system is a two-band metal, as discussed in 
section ||, for t± < 2, and a band insulator for t± > 2. 
Therefore, for J7 = both gaps would be zero for t± < 2, 
and would be set by the band separation 2{t± — 1) for 
t± > 2. For the U = 1 spin gap in Fig. |[ the behavior 
for t± > 2 follows the U — form. However, for t± < 2, 
we expect a nonvanishing spin gap due to the spin liquid 
state, consistent with the exponentially decaying spin- 
spin correlation function. We find numerically that the 
gap in this region is finite, but relatively small, with its 



size associated with the strength of the antiferromagnetic 
correlations in the spin liquid state. We calculate the spin 
gap over a limited range of t± because the numerical ac- 
curacy of the DMRG method and the size of the spin 
gap both become smaller with lower t±. Since one needs 
to calculate small energy differences and do finite size 
scaling in order to accurately determine the spin gap, we 
find that an examination of the ground state spin-spin 
correlations, discussed above, is a more sensitive test for 
the existence of the spin liquid state, given a particular 
numerical accuracy. We omit the charge gap for U = 1 
because the accuracy of the calculated charge gap is com- 
parable to its size in this region. 
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FIG. 3. The spin gap As and the charge gap Ac as a 
function of t\_ at half-filling on a 2 x 32 lattice for various U . 

As U is increased, the relatively sharp transition from 
a spin liquid insulator to a band insulator becomes a 
smooth crossover. This can be understood in terms of the 
spin gap because in the Heisenberg limit, the size of the 
spin gap is set by the perpendicular coupling k, At\/U 
for large J±. One can see the crossover to this behav- 
ior for J7 = 4 and 8, and intermediate However, as 
t± becomes larger, the Heisenberg mapping breaks down 
(when is of the order of U), and the linear growth 
in t± characteristic of the band insulator is restored. For 
large J7, Ac « U/2 at half-filling, and we can see that 
the charge gap is approximately this size for [/ = 4 and 
8. For [7 = 4, the crossover to linear growth of the charge 
gap with is visible for > 2. 

We examine the dependence of the spin gap on U at 
half-filling for the isotropic ladder (i^ = 1-0) in Fig. IJ. 
For small U , the spin gap increases with U as one might 
expect from a weak-coupling picture. For very large C/, 
the spin gap should scale with J k, At'^ /U, since J is the 
only energy scale in the Heisenberg system. Therefore, 
the spin gap will have a peak at some intermediate U, 
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seen to be approximately U = 8 in Fig. |4[ 



< 




U 

FIG. 4. The spin gap As as a function of U at half-filling 
on a 2 X 32 lattice. The inset shows As as a function of filling 
for U = 8. 

Therefore, at half-filling, the two-leg Hubbard ladder 
is an insulator with a spin and charge gap for all U and 
with a crossover from a spin-liquid insulator to a 
band insulator at t±c = 2.0 that grows less sharp as U 
increases. 



B. The Doped System 

The evolution of the spin liquid state of the half-filled 
system with doping is a quite interesting question. As 
discussed above, weak-coupling treatments such as the 
BFRG described above make rather specific predictions 
about the phase diagram of the doped system. In par- 
ticular, it is expected that when the isotropic system is 
doped, the spin gap will remain, there will be one low- 
lying charge mode, and that there will be pairing of the 
charge carriers. In addition, according to the BFRG pic- 
ture, there will a set of phase transitions associated with 
the weak-coupling one-band to two-band transition, and 
with the umklapp processes that occur when one of the 
bands is half-filled. It is not clear to what extent this pic- 
ture will hold within a Hubbard model on a lattice and 
with a realistic band structure, or at the relatively strong 
coupling, U = 8, treated here. The detailed properties of 
the system within each phase, such as the non-universal 
exponents that determine the relative strength of vari- 
ous correlation functions, have not yet been calculated 
in detail using weak-coupling methods and are thus in- 
teresting to determine. 

One can also approach the behavior of this system from 
a strong coupling point of view. This has been done 



within the t-J model by starting fromj-the limit of strong 
magnetic coupling between the chaingl3, J± » J. One 
can then treat the interactions between the rungs per- 
turbatively, and, for small doping, make states that are 
delocalized products of simple excitations of the rungs. 
This gives a physical picture of the behavior for a low den- 
sity of holes, and predicts pairing. In Ref. 13, it is then 
argued that the behavior at isotropic coupling is contin- 
uously connected to that at large large J± by showing 
that, within exact diagonalization calculations, there is 
a continuous evolution of the properties with J±, and 
a qualitative correspondence between the large Jj_ case 
and the isotropic case. A|-«ariational dimer resonating- 
valence-bond (RVB) statec3 used to treat the Heisenberg 
ladder gives a qualitatively similar picture, and helps jus- 
tify the continuity between small and large J±. In the 
following, we will examine the properties of the model 
as a function of filling, and then take slices through the 
phase diagram at constant (n) and varying t± at two 
fiUings: (n) ~ 0.875 and (n) — 0.5 (quarter filling). 

We first examine the behavior of the spin gap and the 
spin-spin correlation function when the half-filled system 
with isotropic coupling {t± — 1.0) is doped with holes. 
Upon doping, the spin gap, as shown in the inset of Fig. 
^ for a 2 X 32 lattice, is reduced but remains finite down 
to a filling of at least (n) — 0.75. Although these results 
are calculated on a finite lattice, we will examine the 
behavior as a function of system size and extrapolate to 
the thermodynamic limit at particular fillings below. In 
the regime shown in the inset of Fig. |4[ the spin gap 
remains finite in the thermodynamic limit. 

The Fourier transform of the q = (q, tt) branch of the 
spin-spin correlation function is shown in Fig. ^ for fill- 
ings between (n) = 0.75 and (n) = 1.0. In order to re- 
duce the effects of the broken translational invariance due 
to the open boundary conditions, we average the spin- 
spin correlation function over up to six i,j pairs for each 
in real space before fourier-transforming. A contin- 
uous range of frequencies is shown because the correlation 
function decayed sufficiently rapidly so that we could ob- 
tain its value for all non-negligible separations. Thus the 
actual lattice size used in the calculation is irrelevant, 
and unlike the usual treatment of a finite system with 
periodic boundary conditions, a straightforward Fourier 
transform yields valid results at any frequency value, pro- 
vided that the correlation function is assumed to vanish 
at larger separations. The same approach could also be 
applied to calculations imposing periodic boundary con- 
ditions, provided that the correlation function decays to a 
value sufficiently close to zero for separations greater than 
half the system size. We show only the q± — n branch 
because the correlations are predominantly antiferromag- 
netic between the chains and the q± = branch is thus 
small and flat, with no interesting features. At half-filling 
((n) = 1.0), the Fourier transform S{q^ tt) peaks at g = tt 
and has Lorentzian line shape, as one would expect for 
a spatially alternating, exponentially decaying function. 
As the system is doped, the peak shifts to smaller val- 
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ues of q commensurate with the lattice filhng, so that it system at t± = 1.0, for which As(0) « 0.12 

,1. 



occurs at the "2A:i^" wave vector 
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FIG. 5. The continuous Fourier transform S{q, vr) of the 
spin-spin correlation function S{i,j,X) for i7 = 8 at various 
fillings on a 2 X 32 lattice. 

We have found that the finite-size effects for the charge 
and spin gaps can be large for a system with open bound- 
ary conditions, even on the largest lattice, 2 x 32 sites. 
We therefore must examine the gaps as a function of sys- 
tem size and extrapolate to the thermodynamic limit. 
We concentrate first on a particular filling, (n) — 0.875, 
relatively close to half-filling. The charge gap Ac and 
the spin gap, are plotted as a function of 1/L for 
three representative t± on lattice sizes up to 2 x 40 in 
Fig. ||. The charge gap. Fig. ^(a), scales approximately 
linearly with 1/L for all three values of t±. The lines 
are least-squares fits to polynomials in 1/L whose or- 
der is chosen so that the number of free parameters is 
one less than the number of L points. The 1/L 
extrapolated value is zero to within the accuracy of the 
extrapolation in all three cases. We believe that this be- 
havior is representative and that the charge gap vanishes 
for all t± and U at this filling. The spin gap. Fig. ^(b), 
shows different scaling behavior for t± = 0.5, 1.0, and 
2.0. For t± — 2.0, the scaling is approximately linear 
in 1/L and the gap goes to zero in the thermodynamic 
limit. For t± = 1.0, the term is substantial, and 

As (1/L ^ 0) « 0.05. For t± = 0.5, the finite-size cor- 
rections are larger, the linear term in 1/L is large, and 
the 1/L^ coefficient is negative. Here As(0) « 0.02, a 
small but finite value. However, the uncertainty is large 
because finite-size corrections are large and because the 
errors in Ag at t± = 0.5 are relatively large. (In cou- 
pled chain systems, the accuracy of DMRG is generally 
inversely related to the coupling between chains.) Also 
shown for comparison is the spin gap for the half-filled 
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FIG. 6. (a) The charge gap Ac and (b) the spin gap As 
for f/ = 8 and (n) = 0.875, as a function of the inverse chain 
length, 1/L. The spin gap at half-filling is shown for compar- 
ison in (b), and the lines are least-squares fits of polynomials 
to the data, as explained in the text. 

The extrapolated 1/L — > values of the spin gap, 
As(0), are shown plotted versus in Fig. |^. For 
tj_ > 1.7, As(0) = 0. For [7 = 0, the one band to two- 
band transition takes place, at ij^ = 1.85 at (ri) — 0.875. 
Within the BFRG pictureE3, the system should be a one- 
band Luttinger liquid, i.e. have no spin gap and charge 
gap for t\_ > t^c- The BFRG also predicts some addi- 
tional phases when tj_ is slightly below the band tran- 
sition point where the Fermi velocities in the two bands 
are very different, and at tj^ = <j^c, when the antibonding 
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band can no longer be treated by linearizing the Fermi 
surface. The additional phases below t^c have gapless 
spin and charge excitations, so the disappearance of the 
spin gap at a somewhat less than tj^c is consistent 
with the existence of these phases, although we cannot 
distinguish between symmetric and antisymmetric spin 
and charge excited states within our DMRG calculation. 
We have, however, not yet seen any evidence for a reap- 
pearance of the spin gap exactly at t±^ — inc- 
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FIG. 7. The infinite system extrapolation of the spin gap, 
As{0), from Fig. 6(b) for (n) = 0.875 and U = 8, and the 
average of the pair-pair correlation function {A{i)A^{j)) av- 
eraged between |i — j| =8 and 12 lattice spacings. 

The infinite system extrapolation of the spin gap, 
As(0), becomes zero at t± ~ 0.38 and negative for 
smaller t±, although is positive for all the finite 
lattices we have examined. In addition, the nearest- 
neighbor spin-spin correlation (M/^Affg) goes to zero 
and has a discontinuity in slope at t± w 0.5, as shown 
in Fig. ||. For U = 0, the bonding band is half-filled 
at t± = 0.4, and the BFRG picture predicts that Umk- 
lapp processes in the bonding band will become relevant 
at this point, leading to a C1S2 phase, consistent with 
a vanishing spin and charge gap. Therefore, the vanish- 
ing of As(0) at t± !=a 0.38 could be associated with this 
weak-coupling feature. 

For tj_ < 0.5, the finite size effects become large and 
the numerical accuracy of the DMRG procedure is re- 
duced, so we have not been able to unambiguously de- 
termine whether or not the spin gap reappears in the 
small phase at this filling, as would be predicted by 
the BFRG. 
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FIG. 8. The nearest-neighbor spin-spin correlation mea- 
sured across a rung, {Mf iMf 2}, at (n) — 0.875 and (n) = 0.5 
as a function of t± for U — 8. Here we have chosen i = 16 on 
a 2 X 32 lattice. 

For the quarter-filled system, the U = Q one-band to 
two-band transition occurs at t± — 1.0. For t± > 1.0, 
the bonding band is half-filled, while for t± < 1.0, both 
bands are fractionally filled. The BFRG predicts that 
within the one-band regime, the system will behave as a 
half-filled Luttinger liquid, i.e. a COSl phase, for which 
Ac ^ and Ag = 0. We exhibit the charge and spin gaps 
as a function of 1/L in Fig. §. For ti_ = 1.2, Ac(0) « 
0.08, whereas As(0) vanishes to within the accuracy of 
the extrapolation. Since the finite size effects are fairly 
large, we have studied lattice size up to 2 x 64. This 
is possible for t± = 1.2 since the numerical accuracy is 
high at this parameter range, with a maximum discarded 
density matrix weight of 6 x 10^^, keeping 350 states on 
the 2 X 64 lattice. For t± = 0.7, which is in the two- 
band region for U — 0, Ac(0) vanishes and As(0) « 
0.07. This behavior is consistent with the CISO phase 
predicted by BFRG in this region. In addition, since 
the bonding band is always less than half-filled for t± < 
t±c at quarter filling, the BFRG calculation predicts no 
additional phases due to umklapp effects for small t±. 
As can be seen in Fig. ||, the behavior of the spin-spin 
correlation across a rung, (Mfj^Mfj), is smooth in the 
small t± region, unlike at (n) — 0.875, consistent with 
the BFRG prediction. 

We now concentrate on the doped spin gap phase 
(CISO), present for isotropic t±, and small to mod- 
erate doping. This phase could be relevant for the 
Lai_a;Sr2;Cu02.5 materijals if, in fact, the ladders are suf- 
ficiently weakly coupledEiJ. The general behavior of holes 
doped into a spin liquid phase might also be relevant to 
the high-Tc materials. In particular, we are interested in 
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whether there is an effective attractive potential between 
the holes, the nature of the state that they form, and the 
type of correlations dominant at long distances. 
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FIG. 9. (a) The charge gap Ac and (b) the spin gap, As, 
plotted as a function of 1/L for f/ = 8, (n) = 0.5, and for 
tx = 1.2 and 0.7. 



We first examine the local behavior as pairs of holes are 
doped into the system. In Fig. |l^ we exhibit the local 
hole density 1 — (n^) plotted as a function of the rung 
index i. Since the system is symmetric with respect to 
exchange of the chains, we calculate the average value on 
a rung. The inhomogeneous structure is due to the open 
boundary conditions as well as the interaction between 
the holes. Due to kinetic energy effects, the boundaries 
repel the holes. When two holes are put into the system, 
as can be seen in Fig. HG, the hole density has one peak 



with a maximum at the center of the ladder. In other 
words, the holes tend to both be near the same rung of 
the ladder, indicating a net attraction between the holes. 
We have calculated the pair binding energy of two holes 
in a half-filled 2 x L lattice, defined as 

Eb = 2So(L, L - 1) - EoiL, L) - E^{L - 1, L - 1) (9) 

on a 2 X 32 lattice with t±_ = 1.0, and obtain a positive 
binding energy, « 0.14. As more holes are added to 
the system, the number of peaks in the hole density is 
equal to the number of hole pairs, indicating that there 
is a tendency for the hole pairs to repel one another. 
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FIG. 10. The local hole density 1 — (jii) as a func- 
tion of position along the chain i for 2, 4, and 8 holes 
((n) = 0.96875,0.9375,0.875) on a 2 x 32 chain for [7 = 8 
and t± — 1.5. 

In order to understand further the nature of the inter- 
action between holes in the antiferromagnetic spin liquid, 
we have examined the pair wave function of two holes di- 
rectly. We do this by calculating the matrix element of 
the pair creation operator, 

{N2\Al,\N^) = {N2\{cl^cl^^-cl^cl,^^)\N,) (10) 

between a state |A^i) with N spin up and N spin down 
electrons, and a state | A^2} with N — 1 spin up and — 1 
spin down electrons. Here the index r — {i, A) denotes the 
position. The pair wave function can then be extracted 
by examining the dependence of this matrix element on 
(r — r'). The result for the matrix element between the 
state with four holes (|A^i}) and the state with two holes 
(|A^2)) at = 1 (TV = 14 on a 2 X 16 lattice) is shown 
in Fig. |ll|. The distance along the ordinate is the rung 
separation plus the cross-chain distance (0 or 1), so that 
a value of 1 corresponds to nearest-neighbor points ei- 
ther along or between the chains, and subsequent points 
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are further along the chains. We have also averaged the 
matrix element for a number of site pairs with the same 
(r — r') in order to reduce the finite-size effects due to 
the open boundaries. 
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FIG. 11. The matrix element {N'i,\(c\ ^c\, ^-c\ ^c\, 
of the pair creation operator between a state |A'^i) with two 
holes and |A^2) with four holes on a 2 x 16 lattice with [7 = 8 
and tx_ = 1.0. The distance i is the rung separation \i ~ j\ for 
intrachain part and \i — i| + 1 for the cross-chain part. 

As can be seen in the plot, the largest amplitude for the 
pair matrix element is for a nearest-neighbor site across a 
rung, and the second largest for a nearest-neighbor site 
on the same chain. The two amplitudes have opposite 
signs, indicating a d-wave-like symmetry. The d-wave- 
like structure of the pair wave function is maintained as 
the separation is increased along the chain. In addition, 
the amplitude is strongly suppressed when the two parti- 
cles are created on sites on the same sublattices, as can be 
seen by the smaller amplitude at every other point. This 
structure can be understood in a strong-coupling picture 
by considering a dimer RVB statecj, which gives a good 
qualitative description of the two-chain Heisenberg sys- 
tem. Here, we have added two holes to a half-filled sys- 
tem that is close to the Heisenberg limit. When two holes 
are added on the same sublattice within the dimer RVB 
state, the RVB state is frustrated in that the number of 
possible valence bond configurations is greatly reduced, 
leading to a higher energy. Therefore, the amplitude of 
the state in which both holes are on the same sublattice 
is suppressed in order to minimize the probability of the 
system being in such a configuration. 

We have also examined the pair wave function by calcu- 
lating similar matrix elements for t_L = 2.0 and — 0.3 
in order to examine the large ty_ Luttinger liquid (ClSl) 
phase and the small t\_ phase. Wc find that the pair 
wave function does not maintain the coherent structure 



shown in Fig. |ll| for separations larger than the nearest- 
neighbor separations, indicating that this structure is 
only present in the CI SO part of the phase diagram. 

We have now ascertained that there is a tendency to- 
wards pairing within the doped, spin gap phase (CISO), 
and have examined the shape of the pair wave function. 
However, in order to determine the nature of the ground 
state, particularly with a view towards possible ordering 
via three-dimensional crossover in the real materials, we 
must examine the long-distance behavior of the corre- 
lation functions. We choose a pairing order parameter 
that creates a pair locally in real space, but has a large 
overlap with the actual pair wave function: 

^^« = 4i.T42i-4ia42T- (11) 

This operator creates a pair in a spin singlet on two neigh- 
boring sites on a rung (labeled i). We then calculate the 
corresponding pair correlation function, (A(i)At(j)), as 
a function of distance between the rungs, |i— i|. This pair 
correlation is shown, plotted on a log-log scale for var- 
ious system sizes and values within the spin-gapped 
phase at (n) = 0.875 in Fig. |l|. 

In order to exhibit the finite-size effects, we show the 
pair correlation function on2x8, 2x16, and 2 x 32 
lattices in Fig. [l^(a). Recall that the local density, as 
seen in Fig. can have substantial dependence on lat- 
tice position. Since the correlation functions on the finite 
lattice with open boundaries are not translationally in- 
variant, we must choose one or more i, j pairs in order 
to best approximate infinite system behavior. Originally, 
we calculated the correlation functions by choosing one 
i, j pair for each |j — j| so that the pair was as sym- 
metrical about the center of the lattice as possible. For 
odd separations, i and j can be chosen symmetrically 
about the lattice center (for a lattice of even length), 
whereas for even separations, they cannot be. Such an 
odd-even effect produced, for example, a spurious peak 
at q = (tt, tt) in the Fourier transform of the spin-spin 
correlation function for the doped system (shown in Fig. 
^). In general, correlation functions in real space calcu- 
lated in this way show spurious oscillations which can be 
seen by comparing the correlation functions at a given 
distance on different lattice sizes. In order to try to re- 
move these effects, we average correlation functions over 
a number of i, j pairs for each \i — j|. We have found 
it best to take enough pairs so that one averages over 
the wavelength of the oscillations in the local quantities. 
Here we average over six pairs, starting with the symmet- 
rically placed pair and then proceeding down the lattice. 
When \i — j| gets near the lattice size, the number of 
possible z, j pairs is smaller due to the proximity to the 
boundaries. As seen in Fig. |l^(a), the pair correlation 
function treated in this way shows only small finite-size 
effects until the last few points before the boundary are 
reached. There are oscillations remaining in the correla- 
tion function, but these are not a finite-size effect and 
are present in the thermodynamic limit. The unphysical 
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peaks at q = (tTjTt) for the doped cases in Fig. ^ dis- 
appeared after this averaging procedure was apphed to 
S{i,j,\). 
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FIG. 12. The pair-pair correlation function {A{i)A^{j)) 
plotted on a log-log scale for for U — 8 and (n) — 0.875. In 
(a) the t± = 1.5 susceptibility is plotted for different system 
sizes, and in (b) the susceptibility is plotted for various t± on 
a 2 X 32 lattice. The dashed line gives the \i — decay of 
the U = pair susceptibility, and the short dashed line in (b) 
shows a decay of \i — ■ 

In Fig. |l^(b), we plot the pair correlation function 
(A(i)A'l'(j)) on a log-log scale as a function of \i — j\ on 
a 2 X 32 lattice for a variety of t± values within the CI SO 
phase at (n) = 0.875. Also shown on the plot are lines 



representing the power law decays 



and 



-Jl 



proximately straight lines on the log-log scale, implying 
a power-law decay of the correlations. If this power law 
has the form \i — jI^"^, v ranges from slightly greater 
than 2 for — 0.8 to approximately 1 for <_l = 1.5. 
In the noninteracting, f7 = 0, system, the pairing cor- 
relation function decays as \i — so the pairing is 
only enhanced over that in the noninteracting system for 
tj_ > 1.0. The strength of the pairing correlations is cor- 
related with the size of the spin gap, as shown in Fig. 
0. The average strength of the pair correlation function, 
calculated by averaging (A(j)A^'(j)) between \i — i\ =8 
and |i— J I = 12 lattice spacings, is shown as a dashed line. 
In the isotropic case, the pairing exponent z/ w 2. Since 
this phase is predicted to be a CISO phase by the BFRG, 
we would expect v = 1/Kp within a Luther-Emery pic- 
ture. Recall that the weakest criterium for dominant 
pair correlations within a Luther-Emery picture is that 
Kp > 1, i.e. that v < 1. Therefore, the pairing corre- 
lations would not be dominant for the isotropic system 
within a Luther-Emery picture, and a direct application 
of this model to the Lai_2;Sr2:Cu02.5 materials would not 
predict a superconducting ground state. 
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FIG. 13. The continuous Fourier transform of the den- 
sity-density correlation function C(q) on a 2 x 32 lattice with 
;7 = 8 and (n) = 0.875 for tj_ = 1.0 and t± = 1.5. The up- 
per curves are the qj_ = n branches and the lower curves the 
q± = branches. 

We now search for a CDW correlation function with a 
reciprocal behavior to the pairing correlations, as in Eq. 
||. In Fig. |l^, we plot the continuous Fourier transform 
of the CDW correlation function 



(12) 



As can be seen, the pair correlation functions are ap- 



at (n) = 0.875, for t± = 1.0 and t± = 1.5, within the 
CISO spin-liquid phase. The real-space correlation func- 
tion is averaged and fourier-transformed in the same way 
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as the spin-spin correlation function of Eq. ^. As can be 
seen, there is no distinct peak at the "2fci?" wave vec- 
tor q = ('7*,7r), where = fcj + fcf. = ■n{n). How- 
ever, there is a peak at q = (7r/4,0) which grows and 
sharpens with increasing t±^ . This wavevector is the cor- 
rect value for the "4fci?" CDW correlations, for which 
|2(fc^ + mod 27r| — 7r/4. Within a bosonization 
treatmentcirciJ, the "4fci?" peak should not appear in the 
density-density correlation function of Eq. |l2[ but only in 
a density-density correlation function composed of four 
fermion densities. We believe that the fact that it ap- 
pears in C(q) is due to the band curvature, which is not 
included in the bosonization picture, and effectively in- 
cludes higher order density-density correlations in C (q) . 

We have isolated the "4^^" portion of the density- 
density correlation function by calculating the correlation 
function of a simple effective rung density 



,cff 



(13) 



Since n° is composed of four fermion operators, rather 
than two as in the ordinary density-density correla- 
tion function, in order to form the correlation function, 
) one must subtract off fourteen different discon- 



cff^cff 



{nT n 



nected pieces, not only {n°^){n°^), as in Eq. 12. We will 



not write all the disconnected terms explicitly here. How- 
ever, the subtraction can be carried out to obtain the real 
space correlation function N{i,j) which can be averaged 
over a number oi i, j pairs and fourier-transformed with 
respect to the rung index, as was done for the spin-spin 
and density-density correlation functions above, to ob- 
tain N{q). The fourier-transformed rung-density corre- 
lation function N{q) is plotted for t± = 1.0 and t± = 1.5 
in Fig. |lj(a). As can be seen, the "4fc/" peak at g = 7r/4 
is now well-defined, and the broad background seen in 
C(q) is no longer present. As in Fig. 13, the peak at 
g = 7r/4 grows in size and sharpness as tj_ is increased 
from 1.0 to 1.5. In Fig. ^(b) we exhibit the averaged 
real-space correlation function N{£), where £ = \i — j\, 
on a log-log scale in order to estimate the exponent of 
the decay. The envelope of the correlations is consistent 
with a straight line on the log-log scale and therefore an 
algebraic decay, although the "4^^" structure makes it 
difficult to directly determine how the envelope decays. 
If we assume a form |« — for the decay of the enve- 
lope, one can estimate 7 to be slightly greater than 2 for 
t± ~ 1.0, and slightly less than 2 for t± = 1.5. 

Therefore, the behavior is opposite to what one would 
expect from the Luther-Emery model: both the pairing 
exponent i' and the "Akp" CDW exponent 7 decrease 
as t± increases within the doped spin-liquid phase, and 
clearly do not obey the relationship v^f = 1. It is not 
clear why the exponents do not behave as expected. One 
possibility is that the identification of this phase with 
the Luther-Emery-like phase found in the bosonization 
calculations is incorrect. Another possibility is that there 
is another CDW-like correlation function that behaves in 
the expected way. A fact that could also be important 



is that the CDW should be quite slowly decaying. For 
example, in the isotropic case, v 2 so that 7 should 
be « 1/2. Perhaps such long-range CDW correlations 
effectively put the system into an ordered CDW state 
on a finite lattice with open boundaries. In this case, 
the "Akp" CDW correlation function, which measures 
the decay of the fluctuations, might not behave in the 
expected way. 
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FIG. 14. The 4fc/ density-density correlation function 
N{i,j), as defined in Eq. (13), (a) fourier-transformed and 
(b) in real space on a 2 x 32 lattice for U = 9>, (n) = 0.875, 
and fx = 1-0, 1.5. 
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III. CONCLUSION 

We have explored the ground state phase diagram of 
the two-leg Hubbard ladder as a function of the band 
filling (n) and the interchain hopping matrix element tj_ , 
concentrating on the case of intermediate to strong on- 
site interaction, primarily U = 8. For the half-filled 
system, we find a spin-gapped insulating phase for all 
t± and U accessible to our numerical density-matrix 
renormalization group calculations. For weak U, there 
is a relatively sharp transition from a spin-liquid in- 
sulator at small t± to a band insulator for t± > 2.0. 
As U is increased, this transition becomes a gradual 
crossover more consistent with a Heisenberg-model pic- 
ture. However, at large enough t±, we would expect a 
return to a band insulator for all finite U. This behav- 
ior is consistent with recent weak-coupling_renormaliza- 
tion group and bosonization calculation£3, as well as, 
in the large U limit, strong-coupling Heisenberg model 
calculationsESlid. 

When the half-filled insulator is doped with holes, we 
find a number of distinct phases. The major feature dom- 
inating the phase diagram is associated with the weak- 
coupling one-band to two-band transition. For t±^ larger 
than the U — band transition value, t±c, our numeri- 
cal results show behavior consist with that of a Luttinger 
liquid, whose band filling corresponds to the filling of the 
U = bonding band. The development of a charge gap, 
expected in the half-filled Luttinger liquid, clearly occurs 
at quarter filling in the two-leg ladder for t± > t±c- For 
t± < ti_c, there is a region near t±c in which the spin gap 
still vanishes. This is consistent with a weak-coupling 
renorrpalization group (RG) calculation of Balents and 
FisheiE3, which predicts phases with gapless spin excita- 
tions where the bands overlap, but the Fermi velocities 
are very different. In addition, for light doping and small 
t^, we find the spin gap vanishes at a finite consis- 
tent with a half-filled bonding band, a point at which 
the Balents-Fisher calculation predicts a phase with a 
vanishing spin gap associated with umklapp processes in 
the half-filled bonding band. For the remainder of the 
t±_ < t±c region, and especially in the light to moderately 
doped isotropic case there is a spin-gapped phase with 
gapless charge excitations in which there are algebraically 
decaying d-wave-like superconducting correlations, and 
"Akp" CDW correlations. However, the pairing correla- 
tions do not seem to be dominant in the isotropic system, 
and exponents of the CDW and pairing correlations do 
not seem to obey theJpverse relationship predicted by a 
bosonization pictureEl'EEl. 
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